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ABSTRACT 

Ramaty et al. (1980) have proposed a model to account for the 5 March 
1979 gamma ray burst in terms of a neutron star corequake and subsequent shock 
heating of the neutron star atmosphere. We elaborate on this model by 
examining the overall energetics and characteristics of these shocks, taking 
into account the e’^-e" pair production behind the shock. The effects of a 
dipole magnetic field in the shock jump conditions are also examined and it is 
concluded that the uneven heating produced by such a field can account for the 
temperature difference between pole and equator implied by the pulsating phase 
of the burst. The overall energetics and distribution of energy between e ~e 
pairs and photons appears to be in agreement with observations if this event 
is at a distance of 55 kpc as implied by its association with the Large 
Magellanic Cloud. 
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I. INTRODUCrr.ON 

If thC5 5 Mjirch 1979 gnmma ray burat is oxtragalactlc as suggested by Its 
association with the supernova remnant N49 in the barge Magellanic Cloud 
(LMG), unique problems are presented in understanding its radiation mechanism 
and energy source. Ramaty ot al. (1980) and Ramaty, hlngenfolter and Bussard 

4 . 

(1981) have proposed an o -o synchrotron cooling aiid annihilation model for 
the radiation mechanism, the primary source of energ>^ being gravitational 
released in the collapse of a neutron star core. This model can naturally 
account for the impulsive phase of the burst, presumably euergiaod by the 
vibrations of the neutron star following the core collapse, v^hich steepen liito 
shocks upon reaching the steep surface density gradient. Such a core collapse 
can easily satisfy the energy requirements oven if the source is in the LMG, a 
view furtljor supported by lt\clependent arguments requiring that the spectrum 
above 500 keV is duo to comptonlzntion of soft photons (Liang 1981).^ 

lu the present paper wo accept the physical association of the event with 
the supernova remnant N49 in the LMG (Roltou 1981) and fvirther expand on the 
model of Ramaty et al. (1980) by studying in greater detail the 
chuu’acterlstlcs of the strong gas shocks presumably responsible for powering 
the impulsive phase of the burst. Assuming complete thermodyi\amlc equilibrium 
behind the shock (including the o’^~e“ pairs), tlte postshock temperatures and 
positron densities are calculated Cor various values of the shock velocity, 
density and magnetic field, appropriate Cor neutron star envelopes. The 
corresponding luminosities and energy distribution between photons and pairs 
are tlion compared to observations and found in reasonable agreement. V'fe 

^Of all Che observed gamma ray bursts, only two have boon identified with 
knovm objects. One is the 5 March event in the LMG and the other, Identified 
recently with Vela X-ray observations (Terrell et al. 1982), Is consistent in 
direction with the binary pulsar SMC X-1 in the Small Magellanic Cloud. 
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furt.hi‘ 1 ' propose? that the subsequent 8 s pulsations are due to the uneven 
surface heating resulting from the neutron star's dipole uvsignotic field. The 
polar regions, having approximately radial fields should experience stronger 
shocks, and Iteucc higher postshock temperatures and lutiinositles than the 
eqviatorlal regions whore the field is approximately parallel to the surface. 
Within the framework of the above assumptions, it is found tlvat reasonable 
values of the field can account for the observed temperature difference 
Inferred from the peaks and valleys of the. pulsating phase. 

An earlier suggestion that quakes within the crust of nearby neutron 
stars power galactic gamma ray bursts was made by Fabian, Icke and Pringle 
(197(i) . This model assumed that quakes, releasing about 10"’^ erg, occur 
within 1 l<m from the neutron star surface and produce galactic, gamma ray 
bursts at a typical distance of 100 pc. 

In Section II the salient features of the burst are critically reviewed 
while In Section III the details and Justification of Che vsssumptlons used in 
tile present paper are given. Section IXIA deals with the question of the 
oorequake ami Che overall energetics associated with it* Xu Section IIIB the 
Rauklne~llugonlot conditions for shocks propagating in the neutron star 
aUiosphore, including the effects of the radiation field, e’*'-e‘“ pairs and 
magnetic field, ace written down and solved numerically. The results are 
summarized in Section IIIC and are used to account for the pulsating phase of 
the burst lu Section IV. Finally in Section V the overall results are 
evaluated and conclusions are drown. 

II, DESCRIPTION OF 5 MARCH GAMMA RAY BURST 
Tlie 5 March 1979 gamma ray burst appears to be of a different class than 
the other gamma ray bursts tlius far reported. The reasons for presuming this 



are the following; 

—A —9 

(1) The peak energy flux in the impulsive phase is •'^.5 x 10 ^ erg cm 
(Mazets et al. 1982), at least an order of magnitude greater than any other 
observed burst. 

(2) Strong positional Identification (Cline et al. 1980, 1982; Felten 

1981) with the supernova remnant, N49, in the LMC places (if this association 

is not coincidental) the source at 55 kpc and implies a luminosity of -1 x 
45 -1 

10 erg s for the initial pulse (Mazets et al. 1982). If the other gamma 
ray bursts are galactic, as implied by the observed frequency distribution 
Rothenflug and Durouchoux 1981, Jennings 1982), the 5 March event is over 4 
orders of magnitude more luminous than the most intense galactic burst and 
produced a total of '■5 x 10^^ erg (Mazets et al. 1979). 

(3) The Impulsive phase of the 5 March event has a very short rise time 
(<2 X 10"'^ s, Cline et al. 1980) and a duration (M).15 s), shorter than most 
other bursts. 

(4) Following the impulsive phase, clear pulses are observed for at least 
three minutes with a period of 8.0 -h 0.05 s. The average pulsed flux is about 
2 orders of magnitude less than the impulsive peak ('•3.6 x 10^^ erg s“^) and 
decreases approximately exponentially with an e-foldlng time of '•50 s (Mazets 
et al. 1979). The total energy emitted during the oscillating phase is '•2 
times as great as emitted in the initial pulse. The pulses show an Interpulse 
of lower intensity that is strongly reminiscent of pulsar pulse profiles and 
implies emission from the magnetic poles of a rotating neutron star. It must 
be pointed out, however, that no observed pulsar has a period as long as 8 s 
and if the 5 March burst is associated with N49, whose age is -10^ years, a 
much shorter period would be expected if an object similar to observed radio 
pulsars produced the burst (Barat et al. 1979). 
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(5) The 5 March event is also unique in that it was followed by three 
other bursts with delays of 0.60, 29 and 50 days (6 March, 4 April and 24 
April) and intensities of 3, 1 and 0.5 percent respectively of the 5 March 
event (Cline 1982). The positional data on these bursts are less precise than 
5 March but are consistent with N49 (more recent observations show a total of 
'•lO bursts from this direction (Mazets 1983)). 

(6) The energy spectrum of the 5 March event is considerably softer than 
typical bursts with most emission at energies below the line feat\ centered 
at about 430 keV . 

In addition to the above unique features, the initial pulse of the- 5 
March event shows a i>road emission line centered at 430 keV with FWHM of -150 
keV (Mazets et al. 1982). This feature has been seen in several of the most 
Intense gamma ray bursts (Mazets et al. 1981) and has been Interpreted as 511 
keV annihilation radiation redshifted in the gravitational field of a neutron 
star with M '•IM^, and R -10^ cm which implies z = (1-2GM/Rc^)“^/^ -- 1 = 0.19 
(Mazets et al. 1979 and Ramaty et al. 1980). The energy flux in the line is 
about 7% of the total impulsive flux (-7 x 10^^ erg s”^, Mazets et al. 1982). 

Ill . THE MODEL 

Fig. FI shows the time history of the 5 March burst. The constrast 
between the very intense, very sharp initial pulse and the slowly decaying 8 s 
pulsations lead us (along with Ramaty et al. 1980) to propose the following 
scenario to account for the main characteristics of this burst. First, the 
energy is released in a phase transition in the neutron star core, releasing a, 

CO 

small fraction of the lO"^ erg in total binding energy as vibrations that 
propagate to the surface at a sizable fraction of the speed of light. Second, 
the initial pulse is generated when the vibrations produce strong, radiation 
dominated, gas shocks in the atmosphere. Third, the slowly decaying pulsating 



7 


OF 

phase is due to the uneven heating of the neutron star surface due to its 
dipole magnetic field. The < 2 x 10“^ s rise time of the initial pulse is 
characteristic of the energy release time of the phase transition. The " 100 
ms decay time of the initial pulse has been shown by Ramaty et al. (1980) to 
be consistent with gravitational radiation damping of non-radlal vibrations in 
a 1 M 9 neutron star. The -^50 s decay time of the pulsations could be either 
the cooling time of the surface layer by conduction to the core and radiation 
at the surface or the damping time of the radial vibrations by the process 
(Langer and Cameron 1969). 

(A) Corequake 

Matter at or near nuclear densities is composed mainly of neutrons. At a 
critical density about twice nuclear density, a new state is thought to exist 
containing a condensed charged pion wave (Hartle, Sawyer, and Scalapino 
1975). This pion condensate has a lower energy per baryon and a significantly 
softer equation of state chan nuclear matter and will affect the mechanical 
properties of neutron stars, such as the mass radius relation. In addition, 
the possibility exists that if the core density slowly increases beyond the 
critical density for pion condensation, a supercompressed, metastable, state 
of neutron matter will be produced (Haensel and Proszynskl 1982 and Haensel 
and Schaeffer 1982) . The slow increase in core density could occur through 
mass accretion,^ by reduced centrifugal forces due to a slowing rotation rate, 
or by cooling (Baym 1981). At some point a phase transition could occur from 


^If accretion produces the increase in core density, Haensel and Schaeffer 
have estimated that more than 10 “^ M 0 must be accreted before the phase 
transition to the pion condensate occurs. If the supernova remnant, N49, is 
'•10^ years old, this implies an accretion rate greater than 10”" Mq yr”^. 

This is over 10^ greater than the upper limit on steady state accretion for 
N49 deduced from X-ray observations (Helfand and Long 1979) and seems to 
exclude accretion as the sole mechanism for increasing the core density unless 
a highly variable accretion process occurred. 
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the uornial (but Bupeteompressed) non-plon condensed state to the pion 
condensed stutot This trmvsition, oi* minlcollnpse, vrould result in a large 
release of energy. 

Haensel and Prosmynski (1982) and Haensel and Schaeffer (1982) have 
calculated in a seini-phouomenological way the basic properties of this 
collapse. These properties depend strongly on the equation of state used and 
due to the uncertainty in models of dense matter are by no means precisely 
known. However, forwhat the above authors consider to be the most realistic 
model, they find the followiugt 

(1) The radius of a collapsing neutron star of -^0.7 decreases by 10 m. 

(2) The upper limit for the time scale of the collapse should be 
approximately the free fall time or (GM/U^) *» 10"^ s for R «=» 10 km, 
consistent with the rapid rise time of the 5 March burst. 

(3) The energy release as estimated by Haensel and Proszyr»ski (1982) 

AO 

would be " 10*^ erg and v>?ould mainly go into heating the neutron star core. 

The proportion tliat v^nould to into vibrational energy in unknovj^n, but if the 
amplitude of the vibration is roughly the change in the radius of the star 
(l.o., '■ 10 m) then the oscillatory energy, can be estimated from 

dimensional considerations (Zel'dovlch and Novikov 1971, p. 366); 

B ^ 10^ (M/MJ (6R/R) erg. (1) 

OSC V o 

*“3 A 7 

If <SR/R " 10 as suggested, the vibrational energy would be 10^^ erg and 

only a small fraction (-^ 0.01) would be sufficient to account for the 

energetics of the burst. For such small amplitudes, the proportion of 

vibrational energy that is lost due to viscosity in the core is insignificant 

and implies tlmt the damping of the vibrations is from the transport of energy 
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to the surface by shock waves, the emission of gravitational waves and 
the I process as mentioned above. 

The supernova remnant, N49, was observipd in soft X-rays 38 days and 2 
years after the burst and only upper limits for the flux were obtained 
(Uelfand and Long 1979, Pizzichini et al. 1982). However, even though 10‘*® 
erg are released by the corequake, it is likely that the only observable 
consequence will result from the small fraction of this energy that is rapidly 
dissipated in the thin surface layer when the vibrations steepen into shocks 
and produces the gamma ray burst. This layer will be thin because, as 
suggested by Fabian et al. (1976), even though the thermal conductivity of 
neutron star material is very high, a shock will not dissipate energy unless 
the shock temperature exceeds the Fermi temperature. This will start to occur 
when the shocks approach the nondegenerate outer crust, i.e., when the density 
drops below - 2.4 x 10® T^^^ gm cm"*® (Lang 1974, p. 253). Even if thermal 
energy is stored at 10^® K (where neutrino losses begin to be important) it 
implies that only densities less than •^2.4 x 10^ gm cm"*® are heated. This 
represents a very thin skin of less than 10^® gm (Soyeur 1980) with an energy 
content much less than the energy of the burst. Therefore, energy must be 
continually supplied to the 'surface by the vibrations and the surface 
temperature will decay with essentially the same decoy time as the 
vibrations. If this is 50 s, the surface will not have been observable in 
soft X-rays 38 days after the burst (Helfand and Long 1979). In addition, 
it is unlikely that the bulk heating of the core would have been observable if 
the event ocurred in the LMC, since 10*® erg distributed among 1 produces 
an increase in particle energy of only nucleon. Tliis implies a 

temperature of less than 6 x 10® K and is at the lower limit of detectability 
in soft X-rays (K. Hurley, private communicaton) only if significant losses 
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did not occur by neutrino emission or gravitational radiation* 

(4) The lifetime of the metastable state depends critically on the core 
density. The Important result is that for cold neutron stars, the transition 
can occur only once and at the point when the core density Increases from 
0.293 to 0.294 fm“^ (Ifm “ 10”^’’^ cm). The age of the star when this occurs 
would depend on how close to the critical density the star was born and on the 
rate at which the core density increased. The phase transition could occur at 
any time from very young to never. In other words, the mlnlcollapse would be 
both very energetic and very rare, a natural explanation for the lack of 
galactic gamma ray bursts of the s.i^e of the 5 March event. 

(5) No disruption of the star results. 

(6) The mlnlcollapse would lower the moment of Inertia by 0.1%, many 

orders of magnitude greater than that implied by radio pulsar glitches, and 
cause a speed-up in the rotation rate (Haensel and Prosz 3 mskl 1982) . The 
subsequent bursts observed at the 5 March location could be the result of 
crustquakes, as envisioned by Fabian et al. (1976), triggered by readjustments 
of the star to its new equilibrium configuration. The total energy released 
in these subsequent burst, though substantial ergs), is a small 

fraction (•^lO”^) of the total burst energy and could conceivably be accounted 
for in this manner. 

(B) Impulsive Phase 

The spectrum observed in the first 4 s of the 5 March gamma ray burst 
(Mazets et al. 1979) and assumed to be indicative of the spectrum of the 
initial 150 ms spike shows a soft, approximately exponential (kT 35 keV) 
spectrum with a clear emission feature, interpreted as redshlfted e'^-e“ 
annihilation radiation, centered at - 430 keV. We propose this emission is 
produced in the strong, radiation dominated, atmospheric shocks that result 
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when the vibrations from the neutron star corequake steepen while propagating 
in the sharp surface density gradient. The shocks will cause the surface to 
expand at a considerable fraction of the speed of light, compressing the gas 
at the surface and producing gas shocks.^ 

(1) Shock jump conditions 

Assuming a reference fram?. where the shock is stationary, the steady 
state jump conditions relating conservation of proton, momentum, energy and 
magnetic fluxes ahead of (subscript 1) and behind (subscript 2) the shock, 
allowing for the production of pairs from the equilibrium photon field behind 
the shock, are as follows: 


>^pl^l = ’"p2“2 


n ,m u, 
pi p 1 








B 


__1 

Sir 


( 2 ) 

(3) 




(protons) 




(electrons) 



+ |aT 


4 

e2 



(positrons) 


(photon and magnetic fields) 


.1 2 
[;rn ,m u, 
'■2 pi p 1 




el 


®x 

“i 


( 4 ) 


^ The speed of sound, s, is determined from s = (dP/dp)^^^^^ where P is the 
pressure. The equation of state for a degenigrat 
is given by (Lang 1974, p. 265); P - 10^ 
g cm“^ we have s - 9 x 10^ cm s~^- 


e, nonrelativlstic neutron gas 
dyn cm“^ and for p 3 x 10^“^ 
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‘rp2”p“2 + ’''”p2‘‘''p2 


1 2 

+ 7j- n -m u„ + Y*n „kT „ 
2 e2 e 2 'e e2 e2 


1 2 

+ n.m u_ + Y'n.kT „ + 2n,m <: 
2+e2 ’e+e2 +e 



h^l “ V2 


(5) 


where 



( 6 ) 


Y 


f 

e 



(7) 


and Y (Yg) is the ratio of specific heats for protons (electrons), y is 
assumed to equal 5/3 and varies from 5/3 to 4/3 depending on the thermal 
kinetic energy of the electrons according to the relation, Yg * i + 
(l/3)[(l+2w')/(l+w')] where w' = mgc2/[ (3/2)kT^2l * 

In the above equations, tip(ng) is the proton (electron) number density, 
n^ is the positron number density, u-^ (U 2 ) is the bulk plasma flow velocity 
ahead of (behind) the shock, ra (m ) is the prot-'n (electron) rest mass, T_ 
(Tg) is the proton (electron) temperature, B is the magnetic field (assumed 
perpendicular to the shock normal), a = 7.56 x 10”^^ erg cm”^ deg”^, and k 
1.38 X 10"^^ erg deg*"^. The terras ^^2 ^gl^^g^l have been omitted 
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from the left hand sides of equations (3) and (4) respectively since 

n^l “ and they are smaller than the corresponding proton terms by the 

factor m„/m„* These terms are included on the right hand sides since the 
c p 

production of pairs allows 0^2 to become greater than Up 2 ** 

The assumptions made in writing tliese equations are the following: (a) 
The shock is locally plane, (b) the gas is ideal and fully ionized, (c) the 
flow velocities, uj^ and U 2 > are nonrelativistic, (d) charge neutrality is 
maintained (i.e. Ug-n^^iip), and (e) equilibrium is established between 
protons, electrons, the photon field, and pairs (i.e. Tp » T^ ® time 

scale short compared to (i) any losses such as radiation (i.e., we assume 
optical thickness » 1 ) and (li) the transit time of the shock through the 
surface density gradient. 

The pair density behind the shock is calculated, assuming equilibrium 
between creation and annihilation at the downstream temperature, Tg 2 > 
integrating over the Ferml~l)irac momentum spectrum (Clayton 1968, p. 274), 

i.e.. 




^2 ^ h ^ exp (Gw + (|)) + 1 


( 8 ) 


n 

where G “ m^c^/(kTg 2 )> w ° the total energy of electrons or positrons in units 
of MgC^ and ff) is the chemical potential in units of kTg 2 * parameter, if), 

is determined iteratively from the conservation of charge ng 2 ~ 114 . “ np 2 
condition, and it is approximately 0 . 1 . 

The above system of equations is solved numerically for given input 
parameters, Upj^, Uj^, Tpj^, and . 

(2) Time scales 

The relevant time scales for this problem are the following: 
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(a) The equilibration time of electron and proton temperatures Is given 
approximately by (Zel'dovich and Ralzer 1966, p. 421) 


O rri3/2 ~1 /QN 

t •“ 8 X 10 8 (9) 

ep 9 27 

27 —3 

where 1127 is the proton number density in units of 10 cm , Tg is the 
electron temperature in units of 10 ^ K and a value of 10 is used for the 
Coulomb logarithm. The equilibration time between electrons and positrons is 
about LOO times less. 

The Compton scattering time is given by 

1 c 1 —I p,”l/2 rj nc 

'• 1.5 X 10 0,J, Tg 2.25 X 10 Tg ' s (10) 

“25 2 

where 0 ^ is the Thomson cross section, ~6.65 x 10 cm . 

(c) The synchrotron loss time (Lang 1974, p. 29) is 


t - 4 X 10"^^b7^ t“^ s (11) 

s 12 9 

12 

where Bj ^2 magnetic field in units of 10 G. 

(d) The annihilation time is 

t - — - 1.3 X 10"^\“1 s (12) 

a oai , 2 / 

+ 

where a is the annihilation rate coefficient and is approximately constant and 

equal to 7.5 X 10“^^ cm^ s“^ (Bussard, Raraaty and Drachman 1979), and n^27 

27 —3 

the positron number density in units of 10 cm 

(e) The plasma frequency is 
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n tn3 1/2 „ 

9 X 10 n Hz 
e 


or 


t ■“ 3.5 X 10 n 8 (13) 

P G<4/ 

27 —3 

where 11^27 is the electron number density in units oC 10 ' cm . 

(f) The electron gyroperiod is 

v"^ - 8 X 10 "^° h~l T 8 a^) 

e 

where Y * (l-(v/c)^) 

(g) The time for the shock to move one gravitational scale height (hg 
600 Tg cm, Liang (1982)), is 

~2 X 10'“ TgJ-l s (15) 

where 3 * u^^/c. This is large compared to the time scales given above and 

justifies the assumption that steady state conditions are obtained as the 
shocks move down the atmospheric density gradient. 

(h) The large temperature gradients produced at the surface will cause 
the shocks to be convective as pointed out by Fabian et al. (1976). They have 
estimated a lower limit on the convective time scale, t^,^, of 

t > 4 X 10“^ s. 
cv — 


( 16 ) 
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The presence of a strong magnetic field is expected to inhibit convention and 
may make t^,^ considerably longer (Proctor and Weiss 1982). Since energy from 
lower hotter layers will be supplied to the rapidly cooling emission layer on 
this time scale, t^^ may be expected to determine intensity fluctuations 
rather than the vibrational frequency or synchrotron cooling time. The fine 
time structure of the 5 March burst has been examined by Barat et al. (1983) 
and some evidence of quasl-periodlc , ^ 25 ras fluctuations is seen. 

The synchrotron loss time is considerably less than the equilibration 
time for electrons and protons implying that the electron temperature 
Immediately behind the shock will be less, due to synchrotron cooling, than 
the proton or photon temperatures contrary to our assumption. However, as 
long as the optical depth remains » 1, any difference in temperature will not 
seriously affect the results of the shock calculation since the shock is 
radiation dominated and, as shown below, the kinetic energy of electrons and 
positrons amounts to less than 5% of the total kinetic energy behind the 
shock. Also, as suggested by Liang (1981), collective effects, whose 
characteristic times, t and v are very short, may be important for 
maintaining equilibrium. In addition, tg << t^ implying, as suggested by 
Ramaty et al. (1980), that pairs will annihilate after losing a significant 
fraction of their kinetic energy by synchrotron emission and hence produce an 
e -e annihilation feature corresponding to a cooler temperature, as observed. 

(3) Input parameters 

g 

(a) At temperatures much below 10 K, the shocks are dominated by the ram 

pressure of the upstream fluid and are essentially independent of the upstream 

9 

temperature. After the atmosphere is heated by the first shocks to 10 K, 
subsequent shocks will be weakened as discussed below. 


(b) We assums^ to be 10^^ G, consistent with pulsar models, and the 
magnetic field configuration to be approximately dipolar. 

(c) The upstream flow velocity, uj^, is essentially the velocity with 
which the neutron star surface expands under the influence of the internal 
vibrations and surface layer shocks. This will depend on the strength of the 
vibrations, the density gradient at the surface, the sound speed in the crust 
material, as well as the gravitational potential and, therefore, the mass and 
radius of the star. None of these quantities are particularly well known, but 
we assume Uj^ to have moderate, nonrelativlstic velocities in the range 4 x 10^ 
“ 10^ km s*"^ considered to be typical for neutron stars. 

(d) The upstream density, is the unshocked gas density at the star 

surface.*^ 

If, for the sake of concreteness, we assume that 1/4 of the surface 
area. A, of a 10 km radius neutron star (A 3 x 10^^ cm^) produces a shock 
with an energy flux, Fgj^, then, for cold upstream material, npj^ and Uj^ must 
satisfy the following condition: 

1 3. „ 

2 pi p 1 sh 

or 

Pj^u^ 6.7 X 10 (cgs units) 

where p, “ n , m . 

1 pi P 


(17) 

U8) 


(C) Results 

Figs. 2, 3 and 4 show the temperature, T^ 2 > positron density, n^, and 


Sjhen the strong Internal shocks reach the surface, the surface layer will be 
vaporized and the atmosphere will consist of whatever gas was present prior to 
the shocks plus the topmost surface layer. 
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compression ratio, r (defined os uj^/u 2 )» plotted against the unshocked 

atmospheric density, Pj^. For example, with Fgj^ -• 10^*^ erg s"^, “ 10^ lun s"*^ 

“-S ?7 -‘3 

and Pj^ ■ 670 g cm , the shock jump conditions yield ■ 9 x 10*' cm T ^,2 " 

1.6 X 10^ K, r ■ 9 and B 2 " 9 x 10^^ G. For the above values, the 

distribution of energy flux behind the shock is divided ns follows (see Table 

1); photons 65%, magnetic field 23.5%, electrons 2.5%, positrons 2%, protons 

2% and the rest mass energy flux of the pairs, 5%. For a given shock energy 

flux, Fgj^, the resulting downstream values of n^ and Tq 2 depend weakly on the 

upstream density, p^, indicating that the shock is radiation dominated. The 

effect of an upstream temperature of 10^ K is also shown in Figs. 3 and 6. 

For a given shock velocity, uj^, the increase in upstream temperature 

results in weaker shocks (i.e., lower compression ratio) and hence leas 

efficient dissipation of the vibrational energy. For upstream temperatures 

greater than the corresponding posts’nock temperatures, T^j^, no shock results 

and the material must cool radiatively before the energy of Che next shock is 

dissipated . 

The effects of the magnetic field on T ^2 c can be significant 
2 2 

when Pj^u^, indicated by the dashed lines of Figs. 2, 3 and A, which 

give the values of the above quantities in the absence of the magnetic 
field. In addition, as mentioned above, the magnetic field is necessary for 
cooling the electrons before they annihilate and hence produce a narrow 
annihilation featvire. 

The observed energy flux emitted in the initial pulse is divided between 
Che e"**-e~ annihilation feature and the continuum and totals about 2 x 10^^ erg 
s"^. the annihilation flux is about 7 x 10^^ erg s~^ yielding a line to 
continuum flux ratio of 0.07. Tlie density of cold positrons needed Co produce 
the annihilation flux is (see equation A) 



19 




l.J 
<1 u 


O 2. , in^3 ~1 

2n,ra c Au« »• 7 x 10 erg s 
T e 2 


( 19 ) 


while Che temperature needed to provide the continuum flux in the photon 
field, assuming blackbody emission is 


A 44 

4aT^2^^2 " ^ (cgs units). 


( 20 ) 


Therefore, the necessary condition to produce the observed line to continuum 
flux ratio is 



( 21 ) 


Table I shows the above ratio for various values of the shock luminosity, 

/ C 1 

Fgj^, and for a reasonable range in shock velocities. For Fgj^ “ 10^^ erg s“'^, 
corresponding to the luminosity of the 5 March event the predicted ratio is 
between 1.3 and 1.5 in good agreement with the value demanded by equation 
( 21 ). 

IV. PULSATING PHASE 

The 8 s pulsations shown in Fig. 1 can be most naturally explained, we 
feel, by emission from the unevenly heated surface (due to the dipole magnetic 
field) of an obliquely aligned, rotating neutron star. The reasons are the 
following: 

(1) The clear pulses with Interpulses precisely 180° out of phase are a 
distinctive signature of emission from a non-uniforraly heated, rotating star 
(Barat et al. 1979, Mazets et al. 1979, Terrell et al. 1980) and indicates 
that global heating of the star has taken place. If the magnetic field axis 
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Is not aligned with the rotation axis, hot magnetic poles will sweep past the 
observer twice for each rotation. If the line of observation is at an angle 
other than 90“ with the rotation axis, one pole will produce a more intense 
pulse than the other. 

(2) The first four main pulses have an average value of kT *= 35.45 keV or 

Q 

T 4 X 10 K (Mazets et al. 1982). If we assume blackbody emission at this 
temperature, an area aligned perpendicular to the line of sight equal to 20% 
of the surface area of a 10 km radius neutron star will have a luminosity of '• 
3,6 X 10^*2 erg s"“^. This is approximately equal to the observed pulsed 
luminosity assuming a distance of 55 kpc. The actual emitting area of the 
star will be somewhat greater depending on the angle between the line of sight 
and a line from the center of the star through the emitting area. 

(3) The difference in temperature between the poles and equatorial 

regions as Inferred from the pulse peak to valley intensity ratio can be 
naturally accounted for by the dipole magnetic field. At the poles, the field 
is approximately perpendicular to the shock surface and no compression of the 
field occurs across the shock. In the equatorial regions the field is 
parallel to the shock surface and will be compressed, absorbing on the order 
kof 50% of the shock energy (see Table 1 for = 10^^ erg s**^, a value 

consistent with the luminosity of the pulsating phase) and resulting in a 
lower post shock temperature. From Fig. 1 we estimate that 

~ 3 - 6 - ( (22) 

eq eq 

implying that 

T 

Y-— ' 0.7 (23) 

poles 
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In Fig. 5 the ratio of dovmstream temperatures with and without a 
magnetic field are calculated for several shock velocities. These 
calculations use the same simplifying assumptions stated above with a constant 

shock luminosity, « 10^^ erg s"^, in equation (18). It is evident that 

the modest temperature difference betweeji the poles and the equator needed to 
account for observations can be easily obtained for field strengths above a 
few 10^^ G. The corresponding shock temperatures are also in the range of the 

Q 

observed values C*- 4.10 K) . 

The problems of energy storage and transfer associated with the impulsive 
phase of the burst are also present in the pulsating phase. The suggestions, 
by Ramaty et al. (1980), of mechanical energy storage and transfer, appears to 
provide, for this phase too, the best solution to the problem. Since 

gravitational radiation must damp all the the non-’radlal modes on short time 

scales (those of the duration of the impulsive phase), the pulsed emission 
would most likely be associated with the energy stored in the radial mode of 
oscillation. The reduced luminosity can then be accounted for in terms of the 
decreased acoustic matching between the neutron star core and its heated 
atmosphere. The acoustic matching is considerably better during the impulsive 
phase, when presumable a large fraction of the oscillatory energy is in the 
higher frequency modes. This energy can be transmitted efficiently through 
the surface density gradient and heat the atmosphere. The transmissivity is 
expected to be significantly reduced for the longer wave length modes of the 
pulsating phase, resulting in less efficient energy transfer and hence lower 
luminosity. 

The magnetic field is expected to play an important role in this 
matching. The radial mechanical motion will shake the magnetic field lines 
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producing outward propagating Alfven waves thus providing a coupling capable 
of tapping the energy of the radial oscillatory modes. Setting aside for the 
moment the question of transmissivity, one can estimate roughly the energy 
flux carried by these waves from F « 5B (B/4ir)cA. Since generally it is 
expected that 6B/B 6R/R, and <SB is not perpendicular to B, the above 
estimate gives 

F “ 1 ^ c 4 xr2 k io‘+ 3 erg 

_3 

assuming 6R/R 10 as argued earlier and using values for B and R typical of 
neutron stars, i.e., 10^^ G and 10^ cm. This value of F is essentially 
identical with that observed in the pulsating phase, provided the source is in 
the LMC. These waves are of course expected to further steepen into shocks 
providing the observed radiation as outlined earlier. 

Eventually the radial oscillations will also damp on much longer time 
scales, due to nuclear processes in the neutron star core (see for example 
Langer and Cameron 1969) thus terminating the burst. Since the energy needed 
in the pulsating phase is provided by shock heating in a way very similar to 
that of the impulsive phase, the emission mechanism for this phase should be 
the same as that for the original spike (Katz 1982). An additional feature of 
the oscillating phase is the difference in temperature between the main pulses 
and the inter pulses. Mazets et al. (1982) report average values of 35.5 and 
31.4 keV respectively. In view of the overall uncertainties of the problem, 
one could possibly attribute this to either asymmetry in the initial 
corequake, favoring heating of a particular area of the neutron star, or 
asymmetry in the star's magnetic field resulting in uneven heating of the 
magnetic poles. (This is also suggested by the pulse profile of Fig. 1, which 



23 


r ■ y 

OF 

shows u deeper mlulmum following the main pulse than following the inter pulse; 
F. C. Michel, private communication). 

V. DISCUSSION AND CONCLUSIONS 

We have examined in greater detail the basic features of the model for 
the 5 March 1979 gamma ray burst suggested by Ramaty et al» (1980), assuming a 
physical association of the burst with the supernova remnant N49. We accept 
the basic features of their model regarding the energy source of the burst 
(corequake of a neutron star), the energy storage (mechanical in vibrations) 
and the transfer of that energy by acoustic waves steepening into shocks in 
the atmosphere. We have in turn examined the characteristics of these shocks 
in greater detail assuming thermodynamic equilibrium behind the shock 
(appropriate for the high density regions below the emission layer) and 
demanding that the energy flux throvigh the shock matches the burst luminosity 
(assumed to be located in N49). 

Assuming values for the density typical to those expected for neutron 
2 4 “3 

star crusts (p =<! 10 -10 g cm ) and velocities of the order of the sound 

4 5 "1 

speed in neutron stars (■^ 4.10 -10 km a ) we obtain shock luminosities which 
are in good agreement with observation, providing that a substantial fraction 
of the neutron star surface radiates. Under the assumption of thermodynamic 
equilibrium we have calculated the corresponding postshock temperatures and 
pair densities for various values of velocity and preshock density. The 
downstream temperature, Tg2> showii in Fig. 2, depends on the velocity, 
density and magnetic field; however T^2 almost independent of these 
quantities for constant F^j^, indicating that the shocks are radiation 
dominated and that most of the shock energy is transferred to the radiation 
field. Therefore, as expected, black body emission at the downstream 
temperature, Tg 2 » produces the observed luminosity if the burst is in the LMC 
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and the radiation lo emitted from a sizat*B fraction of the neutron star 
surfaccv 

The large duty cycle Implied by the time profile of the pulsating phase 
strongly argues in favor of emission over most of the surface of the star and 
in the context of this model suggests that the source is indeed as distant as 
the LMC. If the 5 March 1979 burst luminosity were 10^ times smaller, 
corresponding to a galactic event, involving the whole neutron i3tar surface 
area for the emission as argued above, the corresponding shocks would have 
much lower temperatures and no pairs would be produced. This shock model 
would therefore not be compatible with a nearby event. 

The assumption of thermodynamic equilibrium is of course not valid when 
the shocks break through the atmosphere and the optical depth becomes 1. 
Therefore the present model cannot provide a detailed account of the observed 
spectrum. The detailed interpretation of the spectrum is further complicated 
by the apparent strong temporal evolution during the first 150 ms of the burst 
(Hurley 1983). The model does provide reasonable estimates of the conditions 
and energetics in the denser regions where thermodynamic equilibrium 
prevails. In addition, the ratio of pair annihilation energy to that emitted 
in the continuum is in good agreement with observations for values of 
corresponding to the observed luminosity. If were substantially less, as 
in the pulsating phase or for a galactic burst, no pairs would be produced 
(under the assumptions presently employed). Since no radla^tlve transfer 
effects associated with the last optical depth of the radiation are presently 
taken into account, no strong claims can be made relating these pairs to the 
observed annihilation feature, however the agreement to observations lends 
further credence to the model. 

Finally, this shock heating model provides a natural explanation for the 
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pulsating phase through the enhanced heating of the polar regions due to the 
dipole magnetic field and the rotation of the neutron star. The lower values 
of Inferred from the lower pulsation luminosity again produce temperatures 
consistent with observations if the source is in the LMC. 

The model elaborated upon here, while explaining in general terms the 
properties of this event, involves large scale readjustments in the structure 
of relatively old neutron star cores that must be quite rare, implying that 
some other mechanism powers the more numerous, less energetic galactic gamma 
ray bursts. 
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Fig- I. Initial sharp pulse and first several pulsations of the 5 March 1979 
gamma ray burst. The pulse-lnterpulse structure is evident. Data are from 
Venera 11 and 12 and the figure is taken from Mazets et al. (1982). 

Fig. 2. Downstream plasma temperature, Tg 2 > versus upstream density, pj . 

Solid lines are curves for various shock velocities, uj^ , with = 10^^ G and 

Tpj^ ■=■ Tgj^ = 10^ K. Dashed lines show values of uj^ = 4 x 10^ and 10^ km s~^ 

with = 0 and Tpj^ = T^j^ ■= 10'^ K. The curves labeled (a), (b) and (c) 
correspond to Fgj^ “ 2 10^^, 1 x 10^^ and 5 x 10“^^ erg s“^ respectively. 

Fig. 3. Downstream positron density, n^, versus upstream density, pj^. Solid 
lines are curves for various shock velocities, Uj^, with Bj^ ■=■ 10^^ G and Tpj^ = 
Tgi = 10^ K. Dashed lines show values of Uj^ = 4 x 10^ and 10^ km s"”^ with Bj^ 

- 0 and Tpj^ = T^j^ = 10^ K. The curves labeled (a), (b) and (c) correspond to 
Fgh = 2 X 10^^, 1 X 10^^ and 5 x 10^^ erg s~^ respectively. The dotted line 

shows the effect of increasing the upstream temperature to 10^ K for uj^ ® 10 

km s”^ and Bj^ » 10^^ G. 

Fig. 4. Compression ratio, r, versus upstream density, pj^ . Solid lines are 
for Bjl = 10^^ G. Upper dashed line is for Bj. = 0 and U]^ = 10^ km s“^ while 
lower dashed line is for Bj^ = 0 and Uj^ = 4 x 10^ km s~^ . The curves labeled 
(a) and (b) correspond to ■'sh • 2 X 10^^ and 1 X 10^^ erg s“^ respectively for 
an upstream temperature of 10^ K and B, = 10^^ G. Curves (a') and (b’) are 

Q 

the same as (a) and (b) with an upstream temperature of 10 K. 
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Fig. 5. Ratio of equatorial to polar (equivalent to B “ 0) temperatures 
versus upstream magnetic field, . The shock luminosity, Fgj^ « 1 x 10 erg 
s“^. Dashed line indicates minimum temperature ratio compatible with 
observations. 



TABLE 1 




o 


ON 


o 


VN 

(N 




CM 

O 

X 

oo 


JJ- 

o 

X 

<N 



00 

O 


v^ 

O 

« 

o 


o 

NO 


cs 



T'- 


ON tr> 


rs. 00 

o o 

t • 

o o 


00 VO 
VO VO 


O fvj 
CM 


^ .fM 
O O 


X X 


o o 







S 9E 0 d3d SiNnOO 






(>l) 





COMPRESSION RATIO 




Pihm B 

§P QUALITY 



Oi/1 


ii 


t (Gauss) 











